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ON  THE  BALANCE  EQUATIONS  FOR  A  MIXTURE  OF  GRANULAR  MATERIALS 


S.  L.  Passman 


Introduction. 

The  theory  of  motion  of  a  granular  material  in  three  dimensions  has 
been  given  by  Goodman  and  Cowin  [2],  It  would  be  interesting,  with  certain 
applications  in  mind,  to  extend  this  theory  to  the  case  of  an  arbitrary  finite 
number  of  different  granular  materials  mixed  with  an  arbitrary  finite  number  of 
ordinary  continua  (in  particular,  flu*'ds).  As  a  preliminary  to  this  end,  I  pre¬ 
sent  oalance  laws  for  a  mixture  of  an  arbitrary  finite  number  of  granular  mate¬ 
rials. 

1.  Preliminaries 

The  number  of  symbols  and  the  complexity  of  the  calculations  involved 
in  the  theory  presented  here  is  so  great  that  any  attempt  to  present  either  the 
axiomatic  foundations  or  the  physical  explanations  and  motivations  of  the  theory 
would  enlarge  the  work  greatly.  I  therefore  refer  the  reader  to  the  standard 
v'orks  on  mixture  theory  and  g  anular  materials.  A  list  of  symbols  and  a  short 
set  of  creliminaries  and  standard  equations  is  given  in  order  to  make  the  work 
partial! v  self-contained. 

A  sequence  of  bodies  o  =  is  considered.  A  fixed 

c 

reference  configuration  is  chosen  fot  each  body,  and  in  this  configuration  £ 

a 

is  the  place  occupied  by  a  particle  of  6  .  The  met  ion  of  6  is  the  smooth 

e  a 

mapping 

2$=;s(X,t>,  W:  (1-D 

. .  a  a 

!  in  particular,  I  follow  Truesdell  [1],  and  Goodma:.  and  Cov.-in  [2], 
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of  8  onto  a  region  of  three-dimensional  Euclidean  space  t .  In  general  a 

backward  prime  is  used  to  denote  a  time  derivative  with  #  held  fixed.  Thus 

a 

the  velocity  and  acceleration  of  constituent  a 

h=  at3S  (X  »t)  =£  (25  ,t)  ,  (1.2) 

0  ft  o  a 


25  s  25  (X  » t)  =  x  (js  ,  t)  , 
o  a  a  a 

where  the  second  forms  follow  by  the  assumed  smoothness  of  (1.1).  Often  the 

term  "peculiar _ "  will  be  used  in  place  of  " _ of  constituent  o" 

Thus  &  is  the  peculiar  velocity  of  the  o-th  constituent.  It  is  assumed  that, 
a 

for  each  t  ,  there  is  a  region  of  6  each  point  of  which  is  occupied  simulta¬ 
neously  by  particles  of  each  8 .  Henceforth,  each  formula  written  will  be 
assumed  to  hold  in  subregions  or  at  points  of  this  region. 

Assume  that  each  g  has  a  mass,  which  is  a  measure  on  £  ,  abso¬ 
lutely  continuous  with  respect  to  volume  on  each  configuration  of  8  .  Then  a 

mass  density  p  exists, 
a 

Other  quantities  defined  on  8  are: 

01 

jfc  ,  body  force, 
a 

s  ,  body  heating  , 

T  ,  stress  , 

a 

q  ,  heat  flux  , 

3 

t  ,  internal  energy  , 

k  ,  equilibrated  Inertia  „  (1.4) 


v  ,  volume  distribution 


&  ,  equilibrated  stress  , 

o 

t  ,  equilibrated  body  force, 

g  ,  intrinsic  body  force, 

a 

%  ,  equilibrated  inertial  force  , 

a 

K  ,  inertial  body  force  . 
o 

Growth  terms  are 

c  ,  growth  of  mass  , 

m  ,  growth  of  linear  momentum, 

a 

ifo  ,  growth  of  angular  mom- -n turn  , 

o 

I  X 

e  ,  growth  of  energy  , 

a 

v  ,  growth  of  equilibrated  force, 

£  ,  growth  of  equilibrated  inertia  . 

Here,  lower-case  lightface  symbols  denote  scalars,  lower-case  bold 
face  symbols  denote  vectors,  (with  the  exception  of  jc ),  and  upper-case  bold¬ 
face  symbols  denote  Second -order)  tensors,  thought  of  as  linear  transforma¬ 
tions  on  a  vector  space. 

The  mixture  ma /  also  be  thought  of  as  a  single  body  B  .  Corres¬ 
ponding  to  each  quantity  in  (1,4)  for  8  ,  there  Is  a  similarly  named  and  sym¬ 
bolized  quantity  for  the  composite  body  B  ;  e.  g.  body  force  £  >  body 
heating  s  . 

Truesdeil  [l,3j,  lays  down  three  "metaphysical  principles"  relating 

P  to  B .  Although  it  has  been  ooir.ted  out  that  there  is  some  ambiguity  in  the 

a 
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interpretation  of  these  principles,  I  quote  them,  then  proceed  to  show  that  they 
are  satisfied  by  this  theory.  I  have  some  confidence  that,  for  a  particular 
theory  at  least,  they  may  be  stated  as  axioms  to  be  satisfied  by  the  balance 
laws. 

1.  All  properties  of  the  mixture  must  be  mathematical  consequences 
of  properties  of  the  constituents. 

2.  So  as  to  describe  the  motion  of  a  constituent,  we  may  in  imagi¬ 
nation  isolate  it  from  the  rest  of  the  mixture,  provided  we  allow  properly  for 
the  actions  of  the  other  constituents  upon  it. 

3.  The  motion  of  the  mixture  is  governed  by  the  same  equations  as  is 
a  single  body.^ 

Here  the  "single  body"  considered  in  the  third  principle  is  a  slight 
generalisation  of  the  "granular  material**  cc  ssidered  by  Goodman  and  Cowin 
[*]. 

The  total  mass  density  p  for  the  body  6  is  defined  as  the  sum  of  the 
peculiar  densities 

P=  2  P  ,  (1.6) 

ft 

where  the  symbol  2  ,  here  and  henceforth  is  an  abbreviation  for 

2=5  .  (1.7) 

o  =  1 

The  concentration  g  of  the  o-th  constituent  is  defined  by 

P 


The  principles  are  quoted  from  Truesoell  [1,  p.  83 j,  whore  they  are  subse¬ 
quently  motivated.  It  may  be  noted  that  they  are  somewhat  parallel  in  intent 
to  the  assertion  of  some  molecular  theorists  that  the  properties  of  material 
bodies  are  determined  by  the  p-operties  of  their  molecules. 
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Let  g  be  the  position  vector  from  some  fixed  point  in  £  to  g  . 
Let  <3  ,  &  be  vectors,  let  denote  their  tensor  product,  and  let 


denote  their  outer  product 


(1.9) 


The  densities  of  linear  momentum  and  angular  momentum,  as  well  as 
equilibrated  inertia  and  volume  distribution  momentum  for  B  are  defined  by 


px  =  s  P  *  > 

a  a 


£~px  =  £fi~g 


pk  =  X  p  k  , 
a  3 

pkv  =  2  p  k  v  . 

baa 


(1.10) 


(1.11) 


Here  &  is  interpreted  as  velocity  in  8  . 

2.  Balance  Equations 

I  postulate  the  following  set  of  balance  laws  for  each  constituent  of 


the  mixture,  for  each  part  P  of  0  for  ail  time  t : 

a  a 


/pc  dv  =  (  /  p  d  v)'  , 


(2.1) 


/  P& dv  =  {/  pi  dvr  -  dA  -  /  g  J  dv  >  '2-Z) 

f  ±  2*  &)  =  if  f>  ''PX  dv)'  •  /p^lJQdA  -  /  jg-^pbdv,  (2.  3) 


^There  are  certain  refinements  which  I  eschew  in  order  to  keep  this  work 
relatively  reasonable  in  length.  E.  g.  I  assume  the  existence  of  a  stress 
tensor  rather  than  proving  it  by  the  classical  argument  of  Cauchy. 

#1390  -5- 
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-  $'$•  n  dA  -  .  B  dA 


(s  +  x  •  )2+  v  /  )  dv  , 
0  o  a  a  a 


(2-4) 


/  ptdv  M/pkjdvr  -fft-  B«  -/f  l«  +  9)dv,  (2.6) 

J*p4dv  =  (  J  g  k  av)'  -  ndA-f  p  K  dv.  {2.  o) 

Here  f  dv  denotes  integration  over  the  volume  of  an  arbitrary  part  of  6  and 
J  a 

j>  dA  denotes  integration  over  its  surface  with  unit  outward  normal  n  .  The 
first  four  of  these  equations  respectively  express  growth  of  mass,  linear  momen 
turn,  angular  momentum  and  energy  for  each  component  ot  the  mixture  and, 
except  for  three  terms  in  the  equation  for  energy,  are  of  standard  form.  The 
fifth  equation  expresses  growth  of  equ.librated  force.  An  early  explicit  recogni 
tion  of  a  statement  of  this  nature  as  a  separate  balance  law  is  that  of  Ericksen 
[4 ).  Since  that  recognition,  balance  laws  of  this  type  have  been  widely 
accepted  in  the  theory  of  continua  with  directors.  The  first  explicit  recognition 
of  an  equation  of  the  form  (2.6)  as  a  separate  balance  law  is  apparently  due  to 
Goodman  and  Cowin  [2]  . 

A  standard  argument  applied  to  (2.1)  y ,'iAcS* 

p'csp  +  pdlvj  a  8.p  +  dlv(pj).  (2.7) 

a  a  a  o  1  a  «  o 

the  local  form  of  the  balance  of  mass. 

let  *  bea  scalar,  /actor,  or  tensor  d‘.;i:  -,d  at  points  in  8  ,  By  a 
o  o 

standard  method,  it  may  !;**  shown  that 
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(  j  p  ¥  dv}'  =  j  {  r  $  t  pc  ¥)dv  . 


(2.8) 


If  ^  is  related  to  ^  by 
a 


p*  =  2  p*  , 

aa 


(2.9) 


then  again  by  a  standard  method  it  may  be  shown  that 

p  =  £  p  -  div  £  pijr  a  f  2pc$, 

“  ®  a  o  o  a  o 

where  ii  is  the  diffusion  velocity  of  the  o -th  constituent 
a 

&  =  £-  *  - 
o  a 

By  (2.10),  (1.10)  and  (1.11)  become 

pB  =  ~P*  -  div  Epjssju  tSpcx  , 
a  a  a  a  a  a  a 


(2.10)' 


(2.11) 


p  k  =  2p  k  -  div  £p  k  xj  +  £p  c  k  . 
a  a  a  a  a  a  a 


(2.12) 


pkv-  t  pk  v  r  Sp(k  v  -r  k  v  )  -  div  £p  k  vy  +  £p  ck 


a  ft  a  a  a 


a  a  aa 


Using  the  lemmas  (2.  8)  and  (2.10),  and  noting  that  each  of  (2. 2)-(2.6) 


ho*  for  arbitrary  parts  of  £ ,  I  obtain  the  local  forms 

a 

+  *  -  +  ,,  „ 

PCS  =  p  kS+p&c  -  div  T  -  p  d  , 

a  a  a  aa  aaa 

"t*  x  ^4 

p  (vi  4  {)  Apjj  =  £  Apx^  Apxc-gA  div  J 

a  a  a  a  a  a  a 


(2.13) 


+  1  -  I  -  J0  ^P  b  , 

a  a  a  a 


(2. 1-i) 


Here  "div"  denotes  divergence  with  respect  to  $.  The  divergence  of  a  vector 
is  the  trace  of  its  gradient.  Gradient  is  denoted  by  "grad’',  trace  by  "tr". 
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n 

•  i 


pe  =  p(t+x.x  +  k  v  v  +  ~  k  V2) 


a  a  a  a  a  a  a  a 


a  a 


f  p 


^  ^  ^  <|i 

c  {  i  +  \  x~  +  £  k  v  )  -  2  •  div  X  -  tr(J  grad  $ ) 


a  o 


c  a 


div  a  -  v  div  jb  -  ij  *  grad  v 
o  ft  aft  a 


P  (s  +  x. 
a  a  o 


D 

*■** 

a 


^  )  , 
ao 


(2 .15) 


+ 
p  v 
•  a 


6SST6&£  +  pkjf -d!vfi.g(<tg). 


(2.16) 


X  =pk+pck  -  div  k  -  p  K  . 


a  a  aft 


a  a  a 


(2.17) 


These  local  equations  may  be  reduced  to  simpler  forms  bv  substitution,  yielding 


p  c  =  p  +  p  div  x 
a  a  d  o 


(2.7) 


?h  =  P£  +  P£c  -  div  T  -  p  £  , 

aftft  ftft  ftftft 

+  T 

P  M  =  X  -  I  , 
a  a  o 


+  %  +  i'E,  +  '  +'  +  . »  *  2. 

pe  =  pc  +  pc  (t  -  j  x  )  +  p  m  •  x+pvv  -  p^(j  v  ) 
sac  aft  a  ft  a  aa  aft 


f  p  g  v  -  *r(T  grad  x )  -  div  g  -  jb .  grad  v-ps 
aoa  a  a  aa  Qao 


(2.18) 

{2.  IP) 


4-  (p  K  -  div  k  }  \  'v  Z  . 
a  a  ft  q 


4-  **  %  *  + 

pv  =  pkv  4-pkv  +  pkvc  -  div  h  -  p  (f  -fa)  . 
o  a  a  a  aoa  'ooa  a  aa  c"  ' 


+  v  + 

?  k  spkfpck-  div  k  -  pK  . 
c  ao  a  a  afta 


(2.  201 


(2.21) 


(2.22) 


Equations  (2.7)  and  (2.18M2.22)  appear  to  satisfy  Truesdsll's  second  meta¬ 
physical  principle. 


8- 
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A  assume  that  mass,  linear  momentum,  angular  momentum,  energy, 
equilibrated  force,  and  equilibrated  inertia  are  each  conserved  for  the  mixture, 
that  is 

2  c  =  0  , 
a 

2  h  =  0  , 

a 

2  6  =fl  , 
a 

(2.23) 

2  e  =  0  , 
o  * 

2  v  =  0  , 
a  ’ 

2&  =  0  . 
a 

I  show  that,  through  the  use  of  (2.23),  (2.7)  and  (2.  l3)-(2. 17)  satisfy  Truesdell’ 
third  metaphysical  principle  . 

By  a  standard  argument,  (1.6)  and  (2.7),  (2.  23^  becomes 

p  +  P  div  s  =  0  .  (2.24) 

Define 

T  =  2X  -  2pu®  u  , 
a  o  q  a 

q  =  2g  +  2[JTi}  +  fe  ( v  -  v)  -  p(€  +  iu2  +  kv(v-£v)  y]  , 

a  a  a  a  a  aa  a  aaoo 

h  =  2D  -  2pk  , 

a  a  a  a  a 

=  2k  -  2pkij  , 
a  a  a  a 

pb  =  Sp  b  ,  (2.25) 

o  a 
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ps  =  2ps+2p[b.u+i(v»-v)],  (2.25) 

a  o  a  a  o  a  a  cont  . 

p  t  =  Sp  t  (-  Zj  P  [  u2  +  k  (v  -  v)2] 
a  a  a  a  a  a 

pf  =  2p  i  , 

aa 

pg  =  2P  g  , 
a  a 

pK  =  2p  K  . 
a  i 

These  definitions,  along  with  (1.10)  and  (1.11)  appear  to  satisfy  Truesdell's 
first  metaphysical  principle. 

Summing  (2.18)  over  all  constituents,  noting  (2.23)2  ,  (2.25)^& 
yields,  by  (1.6),  (1. 10)^ ,  and  (2. 12)^, 

pg  =  p  h  +  div  T  .  (2.26) 

Summing  (2. 19)  over  all  constituents,  noting  (2.  23)^  and  (2.25)^  yields 

T  _  jT  (2.27) 

Equations  (2.24),  (2.26),  and  (2.27)  are  respectively  the  usual  local  forms  of 

balance  of  mass,  linear  momentum,  and  moment  of  momentum  for  a  single  body. 

It  is  convenient  to  leave  (2.20)  until  (2.21)  and  (2.22)  are  considered. 

Summing  (2.22)  over  all  constituents  and  noting  (2.23)^  yields 

Spk  +  2p  c  k  -  div  -  2p  K  =  0  .  (2.  28) 

o  a  a  a  a  h 

By  (2.12)  this  becomes 

M 

p  k  -f  div  Spk  u  -  div  £k  ~  2p  K  =  0  .  (2.29) 

a  a  a  a  a  a 

By(2.25)4^1Q  then 

p  k  -  div  k  -  p  K  =  0  .  (2.  30) 


-10- 
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In  the  particular  case  when  k  is  solenoidal  and  K  =  0  ,  (2.  30)  reduces  to 

k  =  0  .  (2.31) 

This  corresponds  to  the  local  equation  for  the  balance  of  equilibrated  inertia 

as  given  by  Goodman  and  Cowin  [2,  equation  (4.9)]. 

Summing  (2. 16)  over  all  constituents  and  noting  (2.  23) _  yields 

5 


2pkv  +  2p£  v+Spkvc  -  div  2h  -  2p  ( I  +  q)  =  0  .  (2.32) 

aN  a  a  a  oaa  a  a  o  a 

By{2.l2)3,  (2.32)  becomes 

pkv  +  pic  v  +  div  2p  k  v  u  -  div  2h  -  2p  (t  +  g)  =  0  ,  (2.  33) 

q  a  a  (T  a  a  a  a 

v/hich,  by{2.25)3  g  ?  is 

pkv  +  pkv  =  div  j}  +  p  [i  +  g)  .  (2.  33) 

This  may  also  be  written  as 

pkH  (div  k+p  K)  v  =  div  &  +  p(f  +  g)  ,  (2.  34) 


where  (2.30)  has  been  used.  In  the  case  where  k  is  solenoidal  and  K  =  0, 

(2.  34)  becomes 

p  k  v  _  div  ij  +  p(i  +  g)  .  (2.  35) 

This  corresponds  to  the  local  equation  for  the  balance  of  equilibrated  force  as 

given  by  Goodman  and  Cowin  [2,  equation  (4. 10)]  . 

The  proof  that  (2.20)  and  (2.23),  yield  the  usual  equation  of  balance 

4 

of  energy,  although  somewhat  repetitive  of  the  derivation  for  the  non-polar  case, 
is  given  in  some  detail  here  in  order  to  emphasize  the  interaction  of  the  new 
terms  in  the  total  internal  energy,  heat  supply,  and  heat  flux  with  the  equi¬ 
librated  stress  and  equilibrated  inertia. 
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Note  that  by  (1. 6),  (1. 10)  and  (2. 11) 


so  tfrat  (2. 43)  becomes 


Sg*  dlvj  +  2tr{TTgradit)  =  div  2TT  u  +  x  *  div  T 
a  a  a  a  o  a 


T  s  v  2, 

+  tr(J  grad  g )  +  div  £p(g®  JS)g  -  div  (p  x  )&  . 


a  a  a 


By  (2.25), 


2  q  =  q  -  S[XTjy  +  &  { v  -  v)  -  p  (t  +  fu2  -  k  v(v  -  \v)  u  ]. 

aaaaaaa 


o  o  a  a  a 

Equations  (1.6),  (1.10)  and  (2.11)  yield 


'2.  . 


.2. 


pu  u  =  Sgx  |  -  22g(x®|c)jc-  2(gx  ) &  +  2px  £ 


Also,  by  (2.  25)?  and  (2.  37) 

2P  «  U  =  2P  «  3  -  2Pt  s 
Qua  a  a  a  a  a 


-  Vj 


‘P  €  x-x(p£  -Sipu2  -  s£pk  (V-  v)2)  , 
a  a  a  a  a  a  a  a 

=  2p«x-p«  i+  ]  Spx2&  -  i(p  x2)x  -  2 1  p  k(v  -  v) 
os  a  a  a 

Expanding  the  iast  term  in  (2.48)  and  using  (1.11)^  2  gives 


aa  a 


Spc  a  =  «  is  -  p  «  £  +  px^  -  £  (p*2)  s 

a  a  a  a  u  a  aa 

.2 


By  (2.  25), 


By  (2.40) 


+  2^pk  v-g  -  ipk  v  x  . 
a  a  a 


2  p  k  wjy  =2)i  -  Ja  • 
a  a  a  a  a 


<  X  =Sp(€  +  ik  ;2)x  -  £{px2)£  -  i(pk  V2)x  . 


‘a '  a  “  a  '  a  a 
Substituting  (2.47)  - (?.  51)  into  (2.  46)  gives 

.T 


SjSl  -  Eft  v+  fo  v  +  £p  *  x  -  P  k  v2  £ 


-5  =  q 

a"  a  o  a  a  aoo 


v  l  ;  pk  l2  #  +  2^(Px2,^-  -rt  (JS®8  )  x 
so  a  a  a  a  a  a  a  o 

.  j  S(Px2)^  . 

0  a 


(2.45) 

(2.46) 

(2.47) 

.  (2.48) 

(2.49) 

(2.  50) 

(2.51) 

(2. 52) 
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By  (2.  25)_  „  ,  (2,25).  becomes 

3,0  O 


2ps=ps-£pb*£+x.pb-2plv+vpf.  (2.  53) 

a  a  a  a  a  aa  a 

The  equation  (2.20)  is  now  reduced  by  summing  over  all  constituents, 

noting  (2.  23)^  ,  and  substituting  values  obtained  from  (2.41),  (2.45),  (2.52) 

and  (2.53).  The  resulting  equation  is 

p  (<  +  k  *  £  +  I  £  »'2  +  k  $  v )  -  tr(TT  grad  x  ) 

-  div  q  -  div  Jtj C-  -  p  (s  f  £  •  b  +  v  t  )  -  0  .  (2.  54) 

Substituting  (2.  26)  and  (2.  32)  into  (2.  54)  yields 
T  2 

p  <  =  tr(J  grad  x  )  +  fc  •  grad  v  +  plcv  +  p  g  iJ  +  div  q  4-  ps.  (2.  55) 
This  is  the  equation  of  balance  of  energy  for  the  mixture.  In  .he  case  when  k 
is  solenoidal  and  K  =  0  ,  by  (2.  31)  it  becomes 

pi  =  tr  (J  grad  jc )  +  •  grad  v  +  p  g  v  +  div  q  +  ps  ,  (2.  56) 

which  agrees  with  the  energy  equation  of  Goodman  and  Cowin  [2,  equation 
(4. 11)].  In  the  case  =  0  ,  it  further  reduces  to  the  classical  equation  for 


the  balance  of  energy. 


3,  The  Entropy  Inequality 


In  addition  to  the  balance  laws  analyzed  in  the  preceding  section,  I 
introduce  a  postulate  of  growth  of  entropy,  analogous  to  the  “second  law  of 
thermodynamics.4'  I  assume,  for  each  constituent  the  existence  of  a  coldness 

$  ,  assumed  to  be  strictly  positive  and  interpreted  as  the  reciprocal  of  the 

a 

absolute  temperature,  an  entropy  g  ,  and  an  entropy  growth  g  .  I  assume  a 

a  a 

balance  law  for  each  constituent  of  the  following  form: 

f  p$l  dv  =  ( f  p  g  tiv)  -  <f  f  q  •  3  dA  -  f  .5  p  s  dv  .  (3.1)^ 

J  q  ^CIG  *  $  ~  'GOG 

By  (2.8)  the  local  form  of  this  equation  is 

4.  v.  4. 

p  g  =pg+pCg-  div  (f  3)  -  *  ps  .  (3.  2) 

0  ao  oa  oqq 

I  introduce  the  axiom  of  dissipation 

2  n  >  0  .  (3.3) 

a 

This  inequality  is  analogous  to  the  conservation  laws  (2.23). 

I  define  the  total  entropy  -g  ,  entropy  flux  6  ,  and  entropy  supply  c- 
mulHplied  by  the  coldness  £  for  the  mixture  by 


p  n  =  2  p  ’i  > 

(3.5) 

G  G 

ii  =  £(.-  q  -  p  g  ij  )  . 

(3.6) 

0  Q  s  5  s 

p^7  =  Sp,?s. 

(3.7) 

1  V  » 

0  q  a 

I  proceed  to  show  that  the  entropy  inequality  for  the  mixture  is  similar  in  form 
to  the  entropy  inequality  for  a  single  continuum  of  the  type  considered. 

^Goodman  and  Cowin  [2]  assume,  for  a  single  constituent,  a  similar  inequality, 
but  with  a  more  general  form  for  the  entropy  “flux"  (the  surface  integral  in  (3.!)). 
They  than  proceed  to  shew  that,  within  a  certain  constitutive  ci-ss.  ;h-.  entropy 
flux  has  the  form  assumed  in  (3. 1).  As  assumption  of  that  nature,  if  made  here, 
would  make  some  of  the  calculations  in  this  section  trivial. 

*It  should  be  noted  that  this  equation  does  not  uniquely  specify  it  <  r  a  . 
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Equation  (3.  5)  is  of  the  form  (2.9).  Therefore  by  (2.10)  , 

P  t)  -  rp  rj  -  d;v  £p  n  jy  +  £p  c  r\  . 
a  o  a  a  a  a  a 


By  ( 3.  2)  this  becomes 


:  j)f,f  aiv  £(p  t]  u  -  a  q  )  -  S  .<  p  s  , 
o  oil  q  s  ~  aaa 


or  ay  (3.6)  and  (3.1)  , 


(3  ;■) 


p  L  r,  :  p  ti  •  QSV  o  -  p.'.  or 
a  ~ 


(3.  *0) 


The  dissipation  axiom  (3.  3)  then  yields 


p  r>  >  div  O  -r  P  A  T  . 


(3. 11) 


This  equation  is  the  same  as  that  of  Goodman  and  Cowin  [2,  equation  (4. 12 ) |. 

It  is  often  convenient  to  state  the  axiom  of  dissipation  in  another  form, 
called  ihe  "reduced  dissipation  inequality,  ”  I  substitute  (2.20)  into  (3.2), 


obtaining 


-  |5+PCTl'(i*  grad  a 


*  f  i  V  2  i  L 

-  *lP«  +  pc(«  '  i-x  )  r  P5»  •  X  -  Pe 
aaa  *  a  a  a  o 

i  '  2  + 

+  p  g  v  *  ;•  k  (p  K  -  div  k  -  p  it) 
oa  ‘aaa  u  a 


-  tr  (I  grad  x  )  -  Jh  •  grad  v  ] 
aaa  « 


Define  the  Helmholtz  free  energy 


so  that 


(3.12) 


(3. 13) 


(3.  14) 
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*i3O0 


Equations  (3.13)  and  (3.14),  when  substituted  into  (3.12),  yield 
+  ,+  v  , .  1^2.-. 
a  a  a  a  a  a  c  a 


+  p  (  —  £  ■  £  v ) 
f  a  a  a 
a 

Tv 

-  jg  •  grad  +  {.  tr(J"  grad  £ )  +  *  fc  •  grad  w 
q  o  ci  q  a  oo  a 


-  {  [  j  v2(p  K  -  div  Jc  -  p  41 )  4  p  §  v]  .  (3- 15) 

g  oao  a  a  a  a  a 

Summing  (3. 15)  ever  all  constituents  and  taking  note  of  (3.  3),  I  obtain  the 

reduced  dissipation  inequality 

£  p.«  [  e  -  m  •  £  -  c  ( 4i  -  \  x2)] 
a  a  a  a  a  a  a 


.  v  ,s 

+  H(-S 

a 


■0 

a  a 


£q 

"a 


grad  $  +  £  <•.  tr(T  grad  £)  +  s  f  J3  '  9rad  v 
ooo  a  Q  o  a 


-  £  [  2  v2(p  K  -  div  };  -  p  p  a  v)]  >  0  .  (3. 16) 

o  aaa  a  a  o  oo 

This  reduces  to  the  result  of  Goodmar  and  Cowin  [2,  equation  (4.15)]  in  the 
case  of  one  constituent. 
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